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ABSTRACT 


Inventory  control  is  an  important  element  of  both  business  and  military  cost  control 
and  readiness.  The  Taiwan  Army  Logistic  Agency  (TALA)  has  used  a  combination  of 
mathematical  inventory  models,  arithmetic  average,  three  month  moving  average  and 
experience  to  project  future  demand.  Implicit  is  that  the  mean  of  monthly  demand  for  an  item 
remains  steady  over  time.  This  assumption  has  proven  to  be  inconect  during  periods  of  force 
reduction,  equipment  is  replaced  or  retired,  or  when  there  is  a  cyclical  demand.  Once  an 
unusual  demand  pattern  occurs,  inventory  control  becomes  impredictable.  Inapplicable 
inventory  methods  in  the  TALA  have  been  estimated  to  cost  as  much  as  several  million 
dollars.  TALA  has  focussed  on  advanced  forecasting  methods.  Exponential  Weighted 
Moving  Average  (EWMA)to  solve  this  problem.  This  may  reduce  inapplicable  inventory  to 
some  extent.  Residual  inventory  and  shortage  are  also  factors  in  cost  control.  In  this  research 
we  will  explore  the  appropriate  approach  to  solve  these  problems  to  make  the  inventory 
control  more  efficient. 
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1.  INTRODUCTION 


A.  BACKGROUND  OVERVIEW 

Taiwan,  Republic  of  China  is  a  small  island  country  located  at  the  western  edge  of  the 
pacific.  The  total  area  of  the  island  is  only  thirty-six  thousand  square  kilometers  (13,971 
square  miles).  The  population  of  this  island  is  more  than  twenty  million  (the  density  is  the 
second  highest  in  the  world).  In  addition,  the  island  has  few  natural  resource  .  Economic 
development  is  very  important  for  Taiwan's  survival  and  prosperity.  However,  the  historical 
threat  from  mainland  China  has  not  decreased  as  predicted.  They  have  maintained  that  they 
will  "never  give  up  resorting  to  forces  to  resolve  Taiwan  issue."  Though  we  have  already 
developed  some  level  of  business  investment  in  mainland  China,  they  still  frequently  hold 
military  drills  on  isolated  islands  near  the  Taiwan  Strait.  So,  the  defense  budgets  still  remain 
at  a  high  level  compared  to  other  developing  coimtries.  This  heavy  burden  will  slow  down 
our  economic  development.  This  is  especially  hard  on  the  democratic  political,  environment 
with  people  eager  for  more  comfortable  life  style.  How  to  maintain  a  balance  of  investment 
between  defense  and  economic  development  has  becoming  a  critical  issue. 

B.  IMPORTANCE  OF  INVENTORY  CONTROL 

There  are  many  different  approaches  that  can  be  used  in  making  tradeoffs  in  the  policy 
of  balancing  defense  and  economic  development;  for  instance,  enhance  training,  reorganize 
the  facilities  or  upgrade  weapon  systems.  But  in  this  study,  I  will  research  the  area  of 
inventory  control. 

Inventory  control  is  a  pivotal  activity  of  any  logistics  organization.  Multi-item 
inventory  systems  encompass  tradeoffs  in  balancing  customer  services  needs  with  operating 
costs.  This  management  task  is  particularly  challenging  in  the  military  setting  where  item 
availability  often  affects  mission  readiness.  Clearly,  one  objective  of  any  inventory  doctrine 
is  success  in  making  decisions  which  minimize  operating  costs  while  providing  an  acceptable 
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level  of  service  for  a  forecasted  rate  of  demand.  The  demand  forecasting  process  and  its 
influence  on  the  inventory  system  effectiveness,  along  with  alternative  inventory  holding  cost 
approximations  are  the  subjects  of  this  study. 

C.  THE  TAIWAN  LOGISTICS  AGENCY  (TALA) 

Before  the  discussion  of  the  various  inventory  control  methodology,  let's  take  an 
overview  of  the  Taiwan  Logistics  Agency  (TALA)  organization  and  its  functions.  The 
organization  of  TALA  is  shown  in  Figure  1. 


Figure  1.  TALA  Organization  and  Function 

1.  Functions  in  the  TALA 
a.  Maintenance  Service 

(1)  Depot  Level.  The  Depot  Maintenance  level  is  the  highest  maintenance 
level  in  the  TALA  organization.  This  level  of  maintenance  performs  the  sophisticated  and 
dedicated  repair  work.  It  performs  inspection,  calibration,  maintenance,  and  repair  of 
equipment  as  well  as  the  training  of  technicians.  Equipment  that  has  been  through  this  level 
of  repair  can  be  regarded  as  equivalent  to  a  new  piece  of  equipment  when  it  is  redistributed 
to  the  units  which  request  it.  Though  this  level's  inventory  control  is  through  BRF  estimated 
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by  skilled  and  experienced  personnel,  a  deviation  from  real  demand  is  unavoidable  due  to 
unexpected  variation  and  human  error.  A  detailed  discussion  will  follow  in  the  successive 
sections. 

(2)  Field  Level.  This  level  is  an  intermediate  maintenance  level.  They 
perform  the  medium  maintenance  and  repair  work.  Most  of  the  work  in  this  level  is  to  change 
the  main  assembly  parts  (engine,  transmission,  differential  gear  box,  wheel  column  and  axle) 
instead  of  repairing  these  assembly  parts  (the  repair  of  assembly  parts  belongs  to  Depot 
Maintenance  Level).  Field  level  also  performs  the  function  of  supporting  the  next  lower 
maintenance  level  (Organization  Level)  with  spare  parts  to  relieve  the  supply  burden  on  the 
lower  level.  Inventory  control  in  this  level  has  less  overhead  compared  to  other  facilities. 

(3)  Organization  Level.  This  level  is  the  lowest  level  of  the  maintenance 
chain.  The  maintenance  work  practiced  in  this  level  is  not  complicated.  The  main  work 
consists  of  lubrication,  adjustment,  cleaning,  and  replacement  of  simple  items,  such  as 
changing  tires,  replacing  light  bulbs,  windshield  blades,  fan  belts,  spark  plugs  and  various 
filters  (air  filter,  fuel  filter  and  oil  filter). 

Though  this  level  only  consists  of  simple  maintenance  work,  the  variation 
of  work  at  this  level  are  much  greater  than  at  other  levels.  Inventory  control  at  this  level  is 
mainly  done  by  less  experienced  personnel  compared  to  other  levels  and  the  estimating 
method  of  spare  parts  is  the  simple  Moving  Average  Approach  (MAA)  to  take  the  variation 
into  consideration.  This  level's  facilities  are  widely  spread  in  the  TALA  system  and  their 
contribution  to  the  total  estimated  forecasting  error  is  considerable.  How  to  improve  the 
MAA  will  be  discussed  in  Chapter  III. 
b.  Supply  Service 

Supply  Services  (SS)  are  composed  of  three  different  functional  groups  described 

as  follows. 

(1)  Inventory  Control.  This  group  is  a  key  functional  area  in  the  SS.  They 
control  the  inventory  level  to  meet  readiness  requirements  and  service  levels,  under  the 
budget  constraint.  They  estimate  the  future  spare  parts  requirements  by  collecting  the 
consumption  data  from  all  over  the  TALA  system.  The  items  that  are  controlled  by  this  level 
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number  more  than  ten  thousand  and  with  average  quantities  in  the  range  of  a  hundred 
thousand.  How  do  the  people  working  here  handle  the  overwhelming  data?  Currently,  they 
use  both  computer  aided  algorithms  (moving  average  approach)  and  experienced  personnel 
to  perform  this  analysis.  This  process  works  well  if  the  average  demand  rate  does  not  vary 
too  much.  However,  in  the  real  world,  there  are  too  many  factors  that  affect  the  demand 
pattern  and  the  accuracy  of  computed  results.  For  instance,  a  periodic  readiness  inspection 
will  make  the  demand  increase  dramatically,  this  will  make  the  moving  average  approach 
less  accurate.  There  is  also  a  potential  human  factor  (fatigue,  sickness,  temper,  overburden) 
that  will  lead  experienced  personnel  to  make  mistakes.  How  to  minimize  these  unavoidable 
mistakes  through  the  forecasting  methodology  is  the  core  of  this  study. 

(2)  Procurement  Group.  The  procurement  unit  working  flow  chart  is  shown 
in  Figure  2.  Basically,  the  procurement  unit  working  process  has  little  relationship  to 
inventory  control  and  other  groups  in  TALA.  The  only  case  for  which  the  procurement  must 
communicate  with  the  inventory  unit  is  when  vendors  offers  a  special  discount  rate  for  a 
different  purchasing  quantity.  This  may  affect  the  inventory  control  process  to  make  the 
inventory  level  higher  or  lower  according  to  the  negotiated  price  and  discoimt  offered. 

(3)  Distribution  Group.  Compared  with  other  facilities  in  TALA,  the 
distribution  group  has  a  simpler  working  procedure.  What  they  are  concerned  with  is 
distributing  spare  parts  at  the  right  quantities;  on  schedule;  to  the  right  place,  and  at  the  same 
time  meeting  the  different  priorities  each  unit  has  set. 
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Request  Time 


Figure  2.  Procurement  Unit  Work  Flow  Chart 


c.  Problems  in  the  Current  TALA  Inventory  Control  System 

As  mentioned  in  the  previous  sections,  we  know  that  out  of  date  computer 
forecasting  approaches  and  the  mistakes  made  by  the  experienced  personnel  working  in  the 
TALA  inventory  system  give  unpredictable  characteristics  to  inventory  system  performance 
and  its  efficiency.  Table  1  shows  the  forecasting  deviation  made  by  the  current  forecasting 
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procedures.  Inadequate  inventory  methods  account  for  five  million  dollars  per  year  or  almost 
fifteen  percent  of  the  annual  inventory  budget. 


Year 

Forecasted  Demand 

Real  Demand 

%  Difference 

0 

14,532 

11,387 

+27.6% 

+1 

13,831 

19,872 

-30.4% 

+2 

45,372 

38,754 

+17.1% 

Table  1.  Forecasting  Deviation  from  Real  Demand 


D.  STOCKAGE 

Stockage  is  the  inventory  level  which  remains  at  the  end  of  each  inventory  cycle.  There 
are  two  different  kinds  of  stockage,  residual  inventory  and  shortage.  Conceptually,  residual 
inventory  is  the  amount  of  material  left  over  or  still  on  hand  when  a  new  order  or  batch 
arrives.  More  formally,  residual  inventory  is  non-zero  when  the  demand  during  leadtime  is 
less  than  the  reorder  point,  in  which  case  it  is  equal  to  the  difference  between  these  two 
values.  At  other  times,  when  demand  during  leadtime  is  less  than  or  equal  to  the  reorder 
point,  residual  inventory  is  defined  as  equal  to  zero.  Shortage  is  the  opposite  of  residual 
inventory.  Conceptually,  shortage  is  the  amount  of  material  back  ordered  because  of 
insufficient  inventory  as  measured  when  a  new  order  arrives.  It  is  also  defined  as  the 
difference  between  demand  during  leadtime  and  the  reorder  point,  but  only  when  the  demand 
value  is  larger;  at  other  times  it  is  zero. 

Superfluous  residual  inventory  will  result  in  waste,  as  represented  by  inventory  holding 
cost.  While  stockouts  may  save  on  holding  cost,  their  negative  impact  on  customer  service 
issues  is  more  profound.  What  are  the  factors  that  affect  stockage  status  and  hence  the 
holding  cost  and  stockouts?  There  are  five  main  factors:  demand  rate,  leadtime,  average 
leadtime,  order  quantity  and  reorder  point.  Two  popular  formulations  for  approximating  the 
holding  cost  area  are  examined  is  this  study  with  respect  to  these  factors  using  the  Crystal 
Ball  spreadsheet  simulation  package. 
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E.  OBJECT  OF  THIS  RESEARCH 


From  the  discussion  in  the  previous  sections,  we  already  know  the  variations  of  various 
factors  will  affect  the  inventory  control.  In  this  study,  several  prevailing  forecasting 
methodologies  will  be  explored  to  determine  which  one  would  minimize  the  total 
forecasting  error  and  eventually  make  the  best  utilization  of  scarce  resources  in  TALA. 

Stockage  is  also  a  key  area  which  will  affect  inventory  control  significantly  and  lead 
to  higher  cost  for  TALA.  Some  of  the  questions  to  be  addressed  in  this  study  are:  to  what 
extent  does  the  underlying  relationships  affect  stockage?  What  are  these  relationships?  How 
to  obtain  the  best  control  over  these  underlying  factors  and  eventually  to  make  the  inventory 
control  more  efficient. 

F.  LIMITATIONS  OF  THIS  STUDY 

As  mentioned  early,  the  number  of  spare  parts  processed  by  TALA  is  greater  than  ten 
thousand.  The  first  limitation  is  that  the  data  retrieved  for  this  study  is  only  a  small  part  of 
the  total  inventory  controlled  by  TALA.  However,  these  data  were  retrieved  randomly. 
According  to  a  normal  distribution  characteristic,  these  sample  data  can  represent  the  major 
characteristic  of  the  whole  TALA  system. 

Second,  the  minor  deviation  factors  might  not  be  detected  by  these  research  due  to  the 
fact  that  unpredictable  variation  factors  cannot  always  be  controlled  as  expected. 

G.  ORGANIZATIONS  OF  THIS  STUDY 

Chapter  II  will  discuss  the  general  forecasting  theorem  and  discuss  the  specific 
function  of  each  theorem.  Chapter  III  will  use  the  data  to  test  each  theorem  discussed  in 
Chapter  II  and  make  comparison  of  each  forecasting  approach  to  specify  which  one  has  the 
potential  merits  over  the  others.  Chapter  IV  will  discuss  the  stockage  analysis  to  draw  the 
underlying  factor's  relationship  in  the  stockage.  Chapter  V  will  present  the  summary, 
conclusions,  and  make  recommendations  as  a  result  of  this  research. 
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II.  FORECAST  METHODOLOGY 


As  stated  in  the  introduction,  forecasting  methodologies  have  two  taxonomies,  the 
qualitative  and  quantitative  techniques.  This  chapter  will  emphasize  quantitative  models  for 
forecasting.  Qualitative  methods  are  used  if  adequate  historical  data  is  not  available  or  cein 
only  be  obtained  at  unreasonable  expense.  That  is  not  the  case  for  the  data  that  is  a  subject 
of  this  thesis.  Supply  usage  data  for  the  TALA  is  readily  available  and  the  collection  methods 
are  in  place  and  have  been  routinely  used  by  TALA  for  several  years. 

This  chapter  is  divided  into  two  parts,  the  first  will  outline  the  forecasting  function;  the 
second  part  will  discuss  time  series  analysis. 

A.  FORECASTING 

A  broad  view  of  forecasting  functions  includes  specific  techniques  and  models  and 
the  significance  of  the  inputs  and  outputs  to/ffom  the  forecasting  model  [Ref.  1].  In  the 
design  of  a  forecasting  function  the  first  step  is  to  specify  its  intended  use  and  the  desired 
output.  The  next  step  is  to  isolate  the  significant  inputs  to  be  used  in  the  model.[Ref.  1] 

1.  Forecasting  Techniques 

Some  facilities  replenish  stock  based  upon  historical  demand.  This  is  reactionary  in 
nature,  and  does  not  normally  reflect  the  timing  of  fiiture  demand.  The  replenishment  system 
purchases  for  a  certain  demand  over  fixed  intervals,  whether  or  not  these  interval  reflects 
actual  demand  over  a  shorter  interval.  This  effectively  supports  the  average  demand  but  the 
greater  the  time  interval  between  when  the  item  is  restocked  and  when  it  is  consiuned,  the 
greater  the  carrying  cost  for  the  inventory.  Furthermore,  if  the  future  demand  is  very  volatile 
or  erratic  in  quantity  and  timing,  the  reorder  quantity  will  be  difficult  to  predict. 

Determining  the  demand  quantity  is  an  important  step  in  ultimately  determining  order 
quantities.  A  common  approach  is  to  simply  take  an  average  for  the  preceding  year.  As  noted 
above,  this  fails  to  consider  the  possibility  of  an  erratic  demand,  seasonal  factors,  and  trends. 
Where  demand  is  computed  monthly  and  future  demand  is  not  expected  to  be  the  same  as 
the  historical  demand,  some  type  of  forecasting  is  required.  It  follows  that  the  better  the 


9 


forecast  of  future  demand,  the  better  that  inventory  levels  can  be  controlled.  Improved 
control  can  translate  into  savings  in  inventory  holding  and  ordering  costs.  Unless  the 
organization  is  willing  to  maintain  excessive  safety  stock  or  incur  frequent  stockouts,  some 
demand  forecasting  is  required. 

There  are  many  different  forecasting  models  that  could  be  used  in  spare  parts  supply 
inventory  management  control.  In  selecting  a  forecasting  method,  we  first  identify  its 
intended  use  and  the  output  results  desired  [Ref  1].  The  next  step  is  to  define  the  inputs 
available  for  the  process,  and  then  the  constraints  and  environmental  factors  affecting  the 
process.  Inputs  could  include  such  items  as  the  demand  history  and  marketing  research 
available.  Constraints  and  environmental  factors  might  include  knowledge  of  special 
situations  which  may  have  affected  the  historical  data,  and  the  availability  of  expert  and 
knowledgeable  personnel.  Outputs  in  the  process  include  the  timing  of  expected  demand, 
broken  down  by  such  segments  as  product,  customer,  and  region. 

2.  Forecasting  Error 

Forecasting  error  is  used  for  two  purposes.  The  first  purpose  is  to  evaluate  the 
effectiveness  of  a  particular  forecasting  method  in  actual  use.  The  second  method  less 
obvious  is  to  make  a  choice  between  forecasting  methods  to  determine  which  method  reflects 
the  least  error  in  selecting  an  organizational  forecasting  method  .[Ref.  1] 

There  are  two  methods  used  to  measure  forecasting  error;  they  include  deviation 
statistics  and  bias  statistics.  Deviation  statistics  give  an  indication  of  the  absolute  magnitude 
of  the  average  error.  Bias  statistics  give  both  direction  and  the  magnitude  of  forec£isting 
errors.  An  example  of  deviation  methods  are  the  mean  absolute  deviation  (MAD)  and  the 
mean  square  error  (MSE).  The  MAD  is  calculated  as  follows: 


MAD 


E  ii;  - 1:1 

jM _ 

n 


(1) 


where: 

=  forecasted  demand  for  period  /, 
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-  actual  demand  for  period  i, 

n  =  number  of  observations  or  time  periods, 

=  deviation  or  forecast  error, 

-  yj  =  absolute  deviation. 

The  MSB  is  calculated  as  follows: 

MSB  ■  — - 

n 


(2) 


Both  MAD  and  MSB  reflect  the  deviation  from  a  forecast  but  the  difference  between  the  two 
methods  is  that  MSB  penalizes  a  forecasting  method  more  for  larger  errors  then  small 
errors.  [Ref  1] 

The  bias  method  of  measuring  forecasting  errors  is  demonstrated  by  the  mean 
error(MB)  and  the  mean  percent  error  (MPB)  methods.  The  MB  is  calculated  as  follows: 


Eci",  ■1’^ 

_M _ 

n 


(3) 


The  MPB  is  calculated  as  follows: 

100  E(y;  -  y)  t  y^ 

MPE  -  - - - 

n 


(4) 


A  positive  bias  indicates  under  forecasting  while  a  negative  result  indicates  that  predictions 
were  too  high  [Ref  1]. 

Statistics  of  forecasting  error  are  required  for  the  following  reasons: 

•  Determine  the  effectiveness  of  a  particular  forecasting  method. 

•  Provide  feedback  for  adjusting  estimates  of  forecast  model  parameters. 
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•  Set  customer  service  level  and  safety  stock. 


3.  The  Forecasting  Process 

The  development  of  a  forecast  should  follow  a  systematic,  six-step  process  which 
includes: 

•  Defining  the  forecasting  problem 

•  Collecting  and  preparing  the  data 

•  Selecting  and  applying  a  forecasting  method 

•  Reviewing  and  adjusting  the  preliminary  forecasts 

•  Tracking  the  forecast's  accuracy 

•  Updating  the  forecast  and  the  forecasting  system 

The  first  step  requires  that  the  forecaster  have  an  understanding  of  the  problem  and  the 
purpose  of  the  forecast.  The  second  step  involves  gathering  the  data  and  ensuring  that  it 
represents  the  information  needed  in  order  to  make  an  accurate  projection.  Step  three 
involves  selecting  the  forecasting  technique  most  appropriate  for  the  available  data,  the 
forecasting  problem,  and  the  situation  or  environment  in  which  the  forecast  is  to  operate. 
Step  four  involves  combining  the  historical  data,  the  forecast  technique,  and  management 
experience  and  judgment  to  evaluate  the  preliminary  forecast  model. 

B.  TIME  SERIES  ANALYSIS 

Time  series  analysis  is  the  forecasting  method  used  in  this  thesis.  It  attempts  to  predict 
future  demand  based  on  past  usage  data.  Time  series  analysis  presupposes  a  relationship 
between  the  observations  of  a  variable  during  successive  and  equal  time  periods.  Time  series 
data  may  contain  up  to  five  components. 

•  Level  Component  or  the  raw  data 

•  Trend  Component  results  from  a  long  term  growth  or  decline 

•  Seasonal  Component  results  from  annual  movements  above  or  below  the  trend 
line 

•  Cyclical  Component  results  from  long  term  cycles,  it  may  or  may  not  be  periodic 
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•  Random  Component  results  from  errors  in  data  collection  and  unusual  natural 
events 

1.  Techniques  for  Time  Series  Analysis 

There  are  several  common  techniques  for  the  analysis  of  time  series  data.  These 
techniques  assume  that  past  data  can  be  used  to  predict  future  data.  Forecasters  should 
understand  that  time  series  analysis  neglects  external  factors  which  may  effect  results  in  the 
future  but  have  not  been  experienced  in  the  past,  and  that  it  will  project  unusual  or  unique 
events  from  the  past  into  the  future.  The  most  common  methods  of  time  series  analysis  are: 

•  Last  Period  Demand 

•  Arithmetic  Average 

•  Moving  Average 

•  Regression  Analysis 

•  Exponentially  Weighted  Moving  Average  (EWMA) 

a.  Last  Period  Demand  (LPD) 

This  method  projects  that  the  next  time  period  demand  equals  the  last  time 
period  so  that  the  forecast  lags  the  actual  demand  by  one  period.  This  method  responds  well 
to  trends  and  is  very  simple.  Some  of  the  disadvantages  are  that  this  method  does  not 
compensate  for  seasonal  and  that  it  will  overreact  to  random  influences. 

Mathematically: 


where: 

=  forecasted  demand  for  period  t, 

,  =  actual  demand  in  the  previous  period. 

b.  Arithmetic  Average 
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This  method  takes  the  average  of  all  past  demand  and  generates  a  projection  of 
future  demand.  This  method  dampens  random  influences  but  does  not  follow  trends  in 
demand  and  does  not  follow  the  seasonal  components. 

Mathematically: 


n 


(6) 


where: 

=  forecasted  demand  for  period  t, 
=  actual  demand  in  period/, 
n  =  number  of  time  periods. 


c.  Moving  Average 

The  method  takes  the  average  of  a  selected  number  of  the  past  demand  data  and 
projects  a  future  demand.  The  selection  is  based  on  the  idea  of  selecting  enough  periods  to 
reduce  the  effects  of  random  influences  and  a  small  enough  number  of  periods  so  that  it  will 
respond  to  the  trend  component. 

The  moving  average  method  is  a  synthesis  of  the  LPD  and  arithmetic  average 
methods.  If  demand  rate  is  smooth  moving  average  works  well  and  responds  quickly  to 
changes  in  demand  but  without  the  fluctuation  of  the  LPD  method.  This  method  does  not 
account  for  the  seasonal  components. 

Mathematically: 


n 


(7) 


where: 
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.  forecasted  demand  for  period  t, 

,  =  actual  demand  in  period  t  -  i, 

n  =  number  of  time  periods  included  in  moving  average. 

d.  Regression  Analysis 

Regression  analysis  establishes  a  relationship  between  various  forecast  variables. 
The  simplest  linear  regression  relationship  is  used  to  fit  a  straight  line  to  the  input  data.  This 
method  compensates  for  trend  effects  but  does  not  follow  seasonal  effects.  Because  of  the 
nature  of  supply  data,  linear  regression  will  not  perform  well  for  those  kinds  of  problems. 
Normal  supply  demand  is  usually  related  in  time  from  one  period  to  the  next.  In  terms  of 
linear  regression  this  would  be  called  autocorrelation.  Autocorrelation  violates  the  conditions 
required  to  produce  a  valid  regression  estimate.  As  an  example,  supply  demand  that  is  very 
high  in  one  period  may  indicate  that  the  next  period  demand  will  be  very  low.  Regression 
analysis  expects  these  kinds  of  fluctuations  to  be  independent  of  each  other. 

e.  Exponentially  Weighted  Moving  Average  (EWMA) 

Exponentially  Weighted  Moving  Average  method  is  a  special  form  of  the 
moving  average  method.  This  method  gives  weight  to  past  data  as  the  moving  average 
method  does,  but  the  weight  varies  geometrically  with  the  increasing  age  of  the  data.  Most 
recent  data  is  weighted  higher  than  older  data.  This  method  is  also  referred  to  as  Exponential 
Smoothing.  The  major  advantage  of  EWMA  is  that  all  historical  data  is  contained  in  the 
previous  projected  demand  and  it  is  not  necessary  to  calculate  on  large  data  collections. 
Mathematically  the  simplest  form  of  EWMA  is: 

X,  .  X,. a{Y,_  ,-i;. ,)  .  a  r,. ,  .  ( 1  -a)  X,_ ,  ,  (8) 


where: 

X,  -  forecasted  demand  for  period  t, 

y,  ,  »  actual  demand  for  period  t-i, 

n  =  Number  of  time  periods  included  in  moving  average. 
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This  is  the  method  used  in  this  thesis  to  predict  supply  demand.  In  addition,  this 
method  has  been  modified  to  incorporate  provisions  for  trend  correction.  To  do  this,  EWMA 
is  used  to  project  the  average  level  and  used  again  to  project  a  level  for  the  trend.  For  this  to 
work  it  is  necessary  to  divide  historical  demand  data  into  two  segments.  Data  in  the  first 
segment  is  used  to  initialize  the  trend  components  in  the  second  segment  of  data  and  linear 
regression  on  the  first  segment  of  data  is  used  to  generate  a  slope  for  the  trend.  The  second 
data  segment  is  used  to  test  various  combinations  of  smoothing  constants.  Once  the 
smoothing  constants  have  been  determined,  future  forecast  can  be  made  using  the  historical 
data  and  parameters  from  the  last  period  in  the  second  data  segment.  EWMA  can  also  be 
used  with  corrections  for  seasonal  components.  For  this  it  is  necessary  to  generate  a  seasonal 
index  that  predicts  demand  level  above  and  below  average. 

If  the  time  series  seasonal  and  trend  components  do  not  exist  the  provisions  for 
those  components  should  not  be  used.  Obviously,  to  answer  this  question,  it  is  necessary  to 
observe  some  historical  data  to  discern  which  components  exist.  This  could  require  the 
storage  of  large  amounts  of  data  which  would  negate  one  of  the  advantages  of  this  method. 
To  avoid  this,  it  is  possible  to  calculate  a  MAD  that  is  modified  for  exponential  smoothing. 
The  MAD  may  be  recalculated  after  each  period  as: 

MAD,  -  d{Y,  -  Y,)  ^  i\-d)  MAD,_^  ,  (9) 


where: 

d  =  exponential  smoothing  constant. 

(1)  EWMA  Extended  for  Trend  Correction.  To  extend  EWMA  to  account 
for  the  trend  components  of  a  time  series  is  necessary  to  maintain  two  adjustments,  average 
level  and  trend. 

Average  level  is  expanded  to  include  previous  trend: 

X,  .  aY,_,  *  (I- a)  i  i;,,  .  r,.,  ).  (10) 
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The  current  forecast  level  is  used  to  calculate  the  trend  as  the  differences 


between  consecutive  forecasts  and  the  trend  component  is  adjusted  as  follows: 


T,  -  biX,-X,_,  )  .  (  1-  i  )  7;.,  . 


(11) 


where: 

Then: 


6  =  the  trend  smoothing  constant. 


T  .  X. 


(12) 


is  used  to  produce  the  forecast  for  the  current  period.  The  EWMA  modified  with  trend 
correction  assumes  level,  trend  and  random  components  are  contained  in  the  time  series. 
For  this  it  is  necessary  to  use  historical  data  to  determine  whether  these  components  are 
present  and  also  to  establish  the  initial  trend  adjustment  r, ,  for  the  model.  [Ref.  1] 

(2)  EWMA  Extended  for  Seasonal  Corrections.  EWMA  reacts  slowly  to 
large  changes  in  demand.  For  normal  services  and  supply  demand  EWMA  can  be  adjusted 
for  a  small  value  of  (a)  and  a  safety  stock  level.  For  a  seasonal  component  it  is  possible  to 
modify  EWMA  to  follow  past  supply  demand  patterns.  The  first  step  is  to  generate  a  set  of 
indexes  that  represent  the  supply  cycle  over  the  length  of  a  season.  These  indexes  will 
represent  demand  above  and  below  average.  The  indexes  are  simply  generated  by  using  the 
total  demand  for  the  season  divided  by  the  number  of  periods.  This  number  the  average 
demand  is  used  to  create  a  number  which  varies  above  and  below  one.  This  number  is 
maintained  as  a  normalized  number  (the  total  is  equal  to  the  number  of  periods).  These 
indexes  are  used  as  an  initial  index  and  represent  the  previous  season.  This  data  can  also 
come  from  an  average  of  several  seasons  in  order  to  reduce  the  effects  of  random 
components  of  the  time  series. 

Actual  demand  data  is  modified  to  remove  seasonal  effects  by  dividing  with 
the  current  seasonal  index.  Mathematically: 
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(13) 


is  multiplied  against  the  index  being  revised  and  all  other  unrevised  indexes  and  will 
maintain  the  seasonal  index  as  normalized  to  the  number  of  periods.  [Ref.  1] 
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III.  COMPARISON  OF  FORECAST  METHODOLOGIES 


A.  INTRODUCTION 

As  discussed  in  Chapter  II,  there  are  various  forecasting  methodologies.  Each  of  them 
have  specific  fimctions  and  can  be  used  to  fit  the  different  demand  patterns:  flat,  trend  and 
seasonal  demand  pattern  or  the  combination  of  the  trend  and  seasonal  demand  pattern. 

In  recent  decades,  a  wide  variety  of  forecasting  methods  have  emerged.  Generally,  they 
can  be  assigned  to  one  of  two  taxonomies:  qualitative  or  quantitative.  Qualitative  techniques 
are  regarded  as  the  more  subjective  of  forecasting  approaches.  Usually  conducted  in  a  setting 
where  historical  data  is  unavailable  or  there  is  little  knowledge  about  forecasting  approaches, 
this  method  often  employs  expert  opinion  in  constructing  a  forecasting  approach.  The 
Taiwan  Army  Logistics  Agency  (TALA)  uses  such  a  method,  called  Best  Replaced  Factor 
(BRF),  when  estimating  initial  stock  levels  for  a  new  item  of  inventory. 

In  contrast,  quantitative  methods  make  extensive  use  of  historical  data.  The  data  serves 
as  an  input  to  various  types  of  mathematics  models  which  compute  the  required  forecast.  Not 
surprisingly,  advances  in  computer  technology  tend  to  popularize  the  quantitative-oriented 
forecasting  methods. 

In  Chapter  I,  we  mentioned  that  most  of  the  time  the  TALA  and  its  Supply  Services 
Division  adopted  a  combination  of  moving  average  and  the  direction  of  experienced 
personnel  to  perform  analysis.  This  is  a  combination  of  qualitative  and  quantitative  forecast 
approaches.  The  qualitative  approach  is  performed  by  experienced  personnel  to  judge  the 
future  demand,  though  it  is  not  as  accurate  as  we  expected,  it  still  is  the  applicable  approach 
when  there  is  no  other  method  available.  The  quantitative  forecast  methodologies  TALA 
uses  most  of  the  time  are  moving  and  arithmetic  average.  The  advanced  forecast 
methodology  of  exponential  moving  weighted  average  (EWMA)  with  trend  or  seasonal 
correction  has  not  been  used  or  researched  by  the  TALA.  The  current  approach  used  by 
TALA  may  be  good  for  some  specific  demand  patterns,  for  example  flat  demand  pattern,  but 
does  not  respond  well  to  the  trend  or  seasonal  demand  patterns. 
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In  this  chapter,  we  will  explore  a  more  advanced  forecast  methodology  and  compare 
it  with  the  methodologies  that  TALA  has  adopted.  We  will  measure  what  benefit  and 
advantages  we  can  obtain  through  EWMA  approach.  The  next  sections  will  discuss  the  flat 
demand  pattern,  trend  demand  pattern  and  the  seasonal  demand  patterns.  In  the  final  section 
will  outline  the  conclusions  from  the  results  of  these  comparisons  and  will  also  discuss  mean 
percent  error  (MPE)  which  will  be  adopted  to  rate  the  accuracy  of  each  methodology. 

B.  FLAT  DEMAND  PATTERN 


Flat  demand  pattern  is  a  demand  with  variation  of  less  than  1 0%  of  the  average 
demand.  Table  2.  shows  a  flat  demand  pattern  collected  from  TALA's  historical  demand 
data .  The  forecast  by  TALA  is  using  the  arithmetic  average  approach. 


Actual  Demand 

Forecasted  Demand  -  TALA's  -  Arithmetic  Average 

115 

130 

110 

123 

108 

118 

109 

116 

113 

114 

111 

114 

114 

114 

113 

114 

117 

114 

117 

114 

108 

114 

114 

114 

MPE -.14% 

MAD-  3.083 

Table  2.  TALA  Flat  Demand  Pattern  Data  Arithmetic  Average 
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We  will  explore  two  other  methodologies,  simple  exponential  smoothing  and  last 
period  forecast  demand(LPD).  LPD  is  a  method  that  simply  assume  that  a  fiiture  period 
demand  is  the  same  as  the  last  period.  The  forecasted  result  is  shown  in  Figure  9  in  Appendix 
A.  and  the  difference  is  shown  in  Table  3. 

We  also  explored  the  simple  exponential  smoothing  approach  in  further  detail  by 
changing  smoothing  constant  a  at  0.01,  0.1,  0.2,  0.3  different  values  trying  to  receive  the 
optimize  smoothing  constant  value  to  make  the  forecast  approach  more  accurate.  The 
comparison  of  each  different  smoothing  constant  is  shown  in  Appendix  A.  Figures  3  and  4 
present  a  comparison  of  arithmetic  average,  LPD  and  EWMA. 


Actual  Demand 

Forecasted  Demand  -  TALA's  -  Last  Period  Forecast 

115 

130 

110 

115 

108 

110 

109 

108 

113 

109 

111 

113 

114 

111 

113 

114 

117 

113 

117 

117 

108 

117 

114 

108 

MPE  =  .01  % 

MAD=  3.25 

Table  3.  TALA  Flat  Demand  Pattern  Data  -  LPD 
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MPE  MAD 


Arithmetic  Average  LPD  Methology 


Figure  3.  Flat  Demand  MAD 


a=  0.1  0.3  0.5  0.8 


Arithmetic  Average  Metholodgy 


Figure  4.  Flat  Demand  MPE 
C.  TREND  DEMAND  PATTERN 

A  trend  demand  pattern  usually  occurs  for  one  or  more  of  the  following  reasons: 

•  Downsize  of  the  military  unit  results  in  decreased  demand. 

•  Deteriorated  or  aged  equipment  results  in  increased  demand. 

•  Subsequently  activated  new  equipment  results  in  a  new  or  increase  of  demand. 
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This  demand  pattern  can  be  detected  by  plotting  the  demand  rate  along  time  intervals. 
If  the  trend  pattern  is  overlooked  or  being  neglected,  the  conventional  forecast  methodology 
(arithmetic  average  or  moving  average)  will  be  wrongly  adopted  by  intuition,  resulting  in 
expensive  forecasted  deviation  from  the  real  demand.  When  we  collected  trend  data  and 
attempt  to  forecast  demand,  using  moving  average  or  arithmetic  average  the  error  from  real 
demand  can  be  approximately  20%. 

An  advanced  alternative  forecasting  methodology  is  exponential  smoothing  average 
with  trend  correction  (EWMA)  is: 

aiY,_,-X,,,)  .  fl  Y,_,  .  {\-a)X,_,  ,  (17) 


where: 


X^  =  forecasted  demand  for  period /. 

y,  ,  =  actual  demand  for  period /-/, 

n  =  Number  of  time  periods  included  in  moving  average. 

With  correction  for  trend  components  : 

X,  =  aY,_,  .  (  1-  a  )  (  i;., .  t;.,  )  ,  (18) 


where: 

T,.,  =  the  trend  from  the  previous  period. 

The  current  forecast  level  is  used  to  calculate  the  trend  as  the  differences  between 
consecutive  forecasts  and  the  trend  component  is  adjusted  as  follows: 

r,  .  6(i;  -  i;.,  )  *  (  1- *  )  T,,  ,  (19) 
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where: 


b  =  the  trend  smoothing  constant. 

Then: 


(20) 


is  the  calculated  forecast  for  EWMA  with  correction  for  the  trend  component  of  supply 
demand  data. 

As  we  explore  this  approach,  we  also  use  the  different  combinations  of  smoothing 
constant  a  and  trend  correction  constant  b  to  optimize  the  forecast  accuracy.  Appendix  B 
present  the  results  for  some  of  the  combinations  of  smoothing  constants  a  and  b.  We  can 
obtain  the  most  accurate  forecast  with  this  approach.  The  deviation  from  real  demand  is 
reduced  from  20%  to  2.84%.  This  is  a  significantly  improved  forecast  accuracy  with  EWMA 
and  trend  correction  to  forecast  demand.  Figures  5  and  6  provide  a  comparison  of  arithmetic 
average,  moving  average  and  EWMA. 

160 
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I  80 
60 
40 
20 

0 
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Figure  5.  Trend  Demand  MAD 
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Figure  6.  Trend  Demand  MPE 

D.  SEASONAL  DEMAND  PATTERN 

As  we  know,  seasonal  demand  will  occur  on  the  big  event  such  as  Christmas  time  in 
the  business  sector.  This  demand  pattern  will  also  occur  in  the  military  sector,  for  example 
during  the  periods  before  cyclic  readiness  inspections  spare  parts  requisition  will  be  higher 
than  normal  times,  and  seasonal  equipment  (anti-slip  chain,  anti-freeze  package),  grade  oil 
(1 OW  for  winter  seasons;  40W  for  summer  season)  demand  will  be  higher  or  lower  according 
to  the  different  seasons. 

How  do  we  forecast  this  kind  of  demand  pattern?  Is  there  any  forecast  methodology 
which  can  forecast  this  kind  of  demand  pattern?  We  explored  the  EWMA  with  seasonal 
correction  approach  to  forecast  the  seasonal  pattern  and  compared  the  result  with  the 
traditional  ones  (moving  average  and  arithmetic  average)  for  TALA  forecast  accuracy. 
Appendix  C  demonstrates  EWMA  adapted  for  seasonal  components  of  TALA  supply  data. 
We  can  see  that  using  the  seasonal  forecast  approach  to  forecast  seasonal  demand  would  be 
better  off  with  a  MAD  of 25.25  while  arithmetic  average  and  a  moving  average  are  at  levels 
of  108.667  and  1 14.083.  Even  a  MAD  of  100  plus  is  misleading  in  that  moving  average  will 
always  lag  the  seasonal  demand  and  arithmetic  average  will  guarantee  that  there  will  be  over 
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MPE  MAD 


and  under  stockage  except  perhaps  two  times  in  a  seasonal  cycle.  Figures  7  and  8  present 
a  comparison  of  arithmetic  average,  moving  average  and  EWMA. 
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Arithmetic  Average  EWMA  Methodology 

Figure  7.  Seasonal  Demand  MAD 
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Figure  8.  Seasonal  Demand  MPE 

E.  CONCLUSIONS 

As  we  know  from  the  previous  sections,  before  we  adopt  a  forecast  methodology  we 
have  to  know  what  the  demand  pattern  is  and  what  time  series  components  it  contains.  This 


can  be  determined  by  plotting  the  historical  demand  rate.  After  deciding  the  demand  pattern 
(flat ,  trend  or  seasonal  demand  pattern)  or  patterns,  then  an  appropriate  forecast  approach 
can  be  adopted  to  create  more  accurate  forecasts  and  reduce  the  deviation  from  real  demand. 
Through  this  approach  we  would  reduce  the  redundant  and  waste  and  reduce  the  defense 
budget  or  allow  expenditures  on  other  requirements. 
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IV  HOLDING  COST 


It  is  important  to  remember  that  holding  cost  is  per  unit  per  unit  time  in  this  model.  For 
any  order  cycle  the  holding  cost  may  be  calculated  by  multiplying  the  area  under  the 
inventory  curve  and  above  zero  by  'h'.  (See  Figure  9) 


I  will  present  two  common  approximations  for  estimating  the  per  unit  time  inventory 
holding  area,  and  hence  holding  cost.  The  differences  between  the  two  are  most  dramatic 
in  the  case  of  order  cycle  with  shortage. 

In  Figure  10,  the  exact  order  cycle  holding  area  is  represented  by  Aiqkl.  The  first 
approximation  is  represented  by  Aijl,  and  the  second  area  approximation  equals  Aiqkl  less 
Akjm.  Clearly  the  first  approximation  would  tend  to  over  estimate  the  actual  value,  while 
the  second  would  under  estimate  the  actual  inventory  holding  area.  In  the  context  of 
inventory  models  such  as  common  EOQ  stochastic  extensions,  algebraic  specifications  of 
the  exact  inventory  holding  area  either  do  not  exist  or  are  prohibitively  cumbersome. 
However,  both  approximations  can  be  readily  specified,  with  the  second  approximation  even 
easier  than  the  first.  The  issue  is  then  to  determine  which  area  approximation  is  more 
accurate,  with  the  realization  that  this  may  depend  on  the  underlying  parameter  values  for 
such  things  as  leadtime  reliability. 
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Figure  10.  Stockout  Approximation  Geometry 


In  both  approaches  we  need  the  expected  value  of  ¥^,  the  average  initial  inventory  of 
the  cycle,  or  any  given  cycle,  the  beginning  inventory  is  equal  to  the  order  quantity  plus  any 
residual  inventory  and  less  any  stockout  from  the  previous  cycle  (assuming  all  stockouts  are 
back  ordered).  Thus  the  average  initial  cycle  inventory  is  equal  to  the  order  quantity  plus 
expected  residual  inventory  and  minus  expected  stockout: 

r,  .  e.r(r)-5(r)  (21) 

where  r  is  the  reorder  point,  and  the  expected  residual  inventory  and  expected  back  order  are 
respectively  defined  by: 

Y(r)  -  ^  (22) 


B(r)  -  /“(*-»•)  fDL(*)  d* 


(23) 


where  fDL(x)  is  the  usage  during  leadtime  probability  density  fimction.  However,  note  that 
the  expected  residual  inventory  minus  the  expected  back  order  is  equal  to  the  reorder  point 
minus  the  expected  demand  during  leadtime  which  in  turn  equals  the  safety  stock: 


Y(r)-B(r)  -  tjopc)  dx  -  r  -  j’xt^(x)  dx 

■  r  -lipL  • 


(24) 
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Therefore,  the  expected  initial  cycle  inventory  is  equal  to  the  order  quantity  plus  safety  stock: 

Q^SS  (25) 


A.  FIRST  APPROXIMATION 

The  first  approximation  involves  specifying  the  average  inventory  level  as  the  simple 
average  of  the  initial  cycle  inventory  and  any  residual  inventory: 

(26) 

With  this  approach,  all  cycles  which  end  in  with  stockout  are  treated  as  thought  they  end  with 
exactly  zero  inventory  --  no  residual  inventory  and  no  back  orders.  The  average  cycle 
inventory  area  is  then: 


1*1.  (“sV  2 


In  this  case  the  per  unit  time  (e.g.  annual)  holding  cost  is: 


HC^‘hY^  =  h  y(r)  * 


This  type  of  holding  cost  approximation  is  adopted  by  Love  [Ref.  2]  and  Brown  [Ref.  3]. 


B.  SECOND  APPROXIMATION 


For  the  second  approximation,  the  average  inventory  level  is  simply  Q/2  plus  safety 


stock: 
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The  somewhat  naive  intuition  behind  this  approach  is  that  inventory  consists  of  two  separate 
parts,  safety  stock  which  is  always  there,  and  the  regular  inventory  which  is  depleted  through 
the  cycle  and  represented  by  Q/2. 

However,  this  approach  ignores  the  fact  that  safety  stock  is  not  always  there,  and  for 
some  cycles  the  corresponding  value  is  negative.  Safety  stock  is  the  average  ending  cycle 
inventory  —  averaged  over  both  residual  inventories  for  some  cycles  and  stockouts  for  other 
cycles.  To  be  more  precise,  safety  stock  is  the  expected  difference  between  the  reorder  point 
and  the  actual  demand  during  leadtime: 

SS  -  fmJix)  dx  (30) 

Normally,  for  most  cycles  this  difference  is  positive,  resulting  in  residual  inventories. 
However,  there  will  always  be  some  cycles  for  which  the  difference  is  negative,  producing 
stockouts.  The  second  approximation  treats  these  back  orders  as  negative  inventory  for  the 
purpose  of  calculating  holding  cost.  That  is,  they  are  treated  as  a  negative  cost,  which  they 
surely  are  not. 

For  cycles  with  stockout,  the  average  cycle  inventory  with  this  approximation  is  an 
average  of  a  positive  number  (Yq)  and  a  negative  number,  the  conditional  expected  stockout: 

*  (31) 


Note  that  the  negative  number  in  this  case  may  be  relatively  large  compared  to  B(r)  because 
the  integral,  representing  the  probability  of  stocking  out,  may  be  quite  small. 

With  the  second  approach,  the  average  cycle  inventory  area  is: 


QY, 


(32) 


and  the  per  unit  time  (e.g.  annual)  holding  cost  is: 
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HC^-hY^^h 


(33) 


This  type  of  holding  cost  approximation  is  adopted  by  several  authors  including  by  Ballou 
[Ref.  4].  The  difference  between  the  two  average  inventory  approximations  is  half  of  the 
expected  stockout: 

r^-T,  -  ^  ^  (34) 

C.  EXACT  INVENTORY  HOLDING  AREA 

The  exact  inventory  area  may  be  determined  by  first  dividing  the  inventory  cycle  into 
two  periods  based  on  when  the  order  is  placed.  Thus  the  first  period  is  the  time  following 
the  arrival  of  an  order  up  until  inventory  gets  down  to  the  reorder  point,  and  the  second 
period  occurs  over  the  length  of  the  leadtime.  The  inventory  area  geometry  of  the  second 
period  will  then  depend  on  whether  the  cycle  ends  with  residual  inventory  or  stockout. 

The  inventory  areas  for  these  periods  are  discussed  below.  An  exact  average  inventory 
area  specification  is  developed  for  the  first  period.  However  for  the  second  period,  the 
specifications  are  not  averages,  but  are  rather  based  on  actual  leadtime  and  demand  values. 
These  second  period  specifications  are  used  through  simulation  to  determine  average  values. 

An  implicit  assumption  of  this  whole  development  is  that  the  demand  in  the  second 
period  is  statistically  independent  of  the  demand  in  the  first.  This  is  an  logical  extension  of 
the  independence  assumption  for  demand  during  any  two  inventory  cycles. 

1.  First  Period  Area 

The  first  period  area  A  is  the  trapezoidal  area  which  is  the  area  between  the  initial 
inventory  level  and  reorder  level,  with  the  base  and  slope  . 
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To  calculate  the  trapezoidal  area  we  can  use  the  variables  specified  above  and 
substitute  them  into  the  trapezoidal  area  equation  (35). 


* 

fr  ' 
A  = 


2 


(Tq  -  r) 
2 


(35) 


So,  we  can  transform  equation  35  into  36: 

^  ^  (Q^SSf  -r^ 

2  Pa 


(36) 


2.  Second  Period  with  Nonzero  Residual  Inventory 

The  second  period  with  nonzero  residual  inventory  is  the  area  demonstrated  in  Figure 
12.  This  figure  represents  the  area  from  the  reorder  point  to  the  end  of  the  cycle  with  a  base 
of  the  leadtime  (LT),  and  with  the  right  hand  side  equal  to  residual  inventory,  i.e.,  the  reorder 
point  minus  demand  during  leadtime  (DL): 
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(37) 


B 


r  + 


-  LT 


X  LT 


where  D  is  the  actual  demand  rate,  LT  is  the  actual  leadtime,  and  DL  is  the  consequent  actual 
demand  during  leadtime  (DL=D*LT). 


Figure  12.  Second  Period 
3.  Second  Period  with  Nonzero  Shortage 

The  second  period  area  with  nonzero  shortage  is  the  area  shown  in  Figure  13  which  is 
after  the  reorder  point  and  before  the  next  order  quantity  has  arrived.  The  inventory  in  this 
area  will  be  depleted  at  some  point  before  the  end  of  the  cycle  which  results  in  a  nonzero 
shortage  during  this  period.  We  can  express  this  area  as: 
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Figure  13.  Second  Period,  Non-zero  Shortage 


C  .  (38) 


But  by  similar  triangles: 


o  LT  J_ 

r  ~  Dl'  D  LT  '  D 


r 


And  hence: 


(39) 


(40) 


D.  SUMMARY 

The  accuracy  of  both  approximations  will  depend  on  several  parameters: 

•  Average  Demand  Rate 

•  Demand  Variability  (Demand  Coefficient  of  Variation) 

•  Average  Leadtime 

•  Leadtime  Variability  (Leadtime  Coefficient  of  Variation) 
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•  Order  Quantity 

•  Reorder  Point 

To  analyze  each  parameter's  effect  on  the  inventory  holding  area  we  utilized  the  crystal 
ball  spreadsheet  simulation  package.  By  changing  each  parameter  value  one  at  a  time,  and 
keeping  other  parameters  unchanged,  we  can  observe  each  parameter's  effect  on  the  holding 
area  for  both  the  first  and  second  approximations,  as  well  as  the  exact  formulation.  For  all 
simulations,  both  leadtime  and  demand  were  assimied  to  have  log-normal  distributions.  Each 
simulation  included  1 ,200  trials. 

Initially  all  parameter  values  were  specified  as  in  Table  4.  Average  demand  was 
always  kept  at  this  value  (1,000  imits  per  unit  time).  However,  all  of  the  other  parameters 
were  varied  to  determine  if  the  relative  accuracy  of  the  two  approximations  depended  on 
these  parameter  values.  First,  demand  variability,  leadtime  variability  and  average  leadtime 
were  independently  varied  as  indicated  in  Table  5.  Then,  order  quantity  and  reorder  point 
were  varied  together  as  presented  in  Table  6.  For  example,  the  bottom  row  indicates  that  a 
reorder  point  value  of  5,000  units  was  examined  in  conjunction  with  two  order  quantity 
values,  5,000  and  10,000. 

Figure  14  presents  the  results  of  this  analysis  for  the  case  where  demand  variability  was 
varied.  These  results  are  typical  in  the  following  respects: 

•  As  expected,  the  simulated  exact  average  area  was  always  less  than  the  first 
(Love/Brown)  approximation. 

•  Contrary  to  expectations,  the  simulated  exact  average  area  was  sometimes  less 
than  and  sometimes  greater  than  the  second  (Ballou)  approximation 

•  In  absolute  terms  (ignoring  signs),  the  second  (Ballou)  approximation  was  always 
superior  to  the  first  (Love/Brown)  approximation,  usually  by  an  order  of 
magnitude  or  more.  Occasionally  they  provided  approximately  the  same  accuracy. 

The  results  for  the  other  variations  (2nd  and  3rd  columns  of  Table  5,  and  Table  6)  are 

presented  in  Appendix  D. 

These  results,  as  illustrated  by  Figure  14,  merit  further  investigation.  In  particular,  the 
imexpected  relationship  between  the  second  (Ballou)  approximation  and  the  simulated  exact 
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values  needs  to  be  explained.  The  next  logical  step  would  involve  using  the  real  analysis 
(calculus)  associated  with  probability  theory  to  investigate  these  questions. 


Average  Demand 

1000 

Demand  CoefRence  of  Variation 

0.01 

Average  Leadtime 

1 

Leadtime  Coeffience  of  Variation 

0.01 

Reorder  Point  (ROP) 

1000 

Qrder  Quantity  (Q) 

5000 

Table  4.  Primary  Parameter  Values 


Demand  Coefficient  of  Variation 

Leadtime  Coefficient  of  Variation 

Average  Leadtime 

0.01 

0.01 

0.001 

0.05 

0.05 

0.01 

0.1 

0.1 

0.1 

0.2 

0.2 

0.2 

0.3 

0.3 

0.5 

0.5 

0.5 

0.75 

0.8 

0.8 

1.0 

1.0 

1.0 

1.1 

1.2 

1.2 

1.5 

1.5 

2.0 

2.0 

5.0 

5.0 

10.0 

10.0 

Table  5.  Leadtime  and  Demand  Parameter  Values 
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Difference  of  Area(%) 


Reorder  Point  (ROP ) 

Order  Quantity(Q) 

1 

1 

10 

100 

1000 

5000 

10000 

100 

100 

1000 

5000 

10000 

1000 

1000 

5000 

10000 

5000 

5000 

10000 

Table  6.  Combination  of  ROP  &  Q 


Figure  14.  Demand  Coefficient  of  Variation  V.S.  Inventory  Area 
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V.  SUMMARY,  CONCLUSIONS  AND  RECOMMENDATION 


A.  SUMMARY 

The  TALA'S  inventory  management  system  has  been  historically  ineffective  in 
preventing  the  accumulation  of  inapplicable  inventories  when  equipment  is  retired  from 
service  and  has  failed  to  anticipate  shortage  due  to  increasing  demand  when  new  equipment 
has  entered  service.  Both  of  these  shortcomings  are  due  in  large  part  to  inappropriate  or 
nonprofessional  forecasting  approach.  The  moving  average  or  arithmetic  average  forecasting 
methodology  will  only  perform  well  for  some  special  demand  patterns  (in  most  cases  the  flat 
demand  pattern).  The  flat  demand  pattern  rarely  occurs  and  these  forecasting  approaches 
caimot  provide  reliable  demand  prediction.  There  are  several  advanced  forecasting 
methodologies  designed  for  tracking  complicated  demand  patterns  (trend  demand  pattern, 
seasonal  demand  pattern  and  the  combination  of  trend  and  seasonal  demand  pattern).  These 
Exponential  Weighted  Moving  Average  (EWMA)  methodologies  [Ref.  1.  pp.  53-63]  can  be 
applied  efficiently  if  the  demand  pattern  for  the  item  has  been  researched  by  plotting  the 
items  historical  demand  rate  along  a  time  line.  Once  the  demand  pattern  has  been 
established,  we  can  use  the  relevant  forecast  methodology  (simple  EWMA,  EWMA  with 
trend  correction,  EWMA  with  seasonal  correction)  to  project  the  future  demand. 

In  this  research  we  collected  some  TALA's  historical  forecast  demand  data,  analyzed 
those  items'  demand  patterns,  and  reforecasted  them  at  a  projected  demand  based  on 
EWMA.  The  results  shows  there  was  an  improved  forecast  accuracy  using  EWMA.  If  TALA 
would  adopt  this  approach,  they  may  save  millions  of  dollars  per  year  and  improve  the 
readiness  and  efficiency  of  their  inventory  system. 

As  mentioned,  there  are  two  different  kinds  of  stockage,  residual  inventory  and 
shortage.  Superfluous  residual  inventory  will  result  in  waste,  as  represented  by  inventory 
holding  cost,  and  while  stockout  may  save  on  holding  cost,  there  is  a  negative  impact  on 
customer  service  and  readiness.  There  are  complicated  factors  that  interact  and  effect  the 
stockage  level  and  inventory  control.  The  most  influential  are  average  leadtime,  leadtime 
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variation,  average  demand  rate,  demand  variation,  reorder  point  and  order  quantity.  Two 
popular  formulations  which  take  these  factors  into  consideration  to  compute  the  holding  cost 
area  were  examined.  The  first  approximation  involves  specifying  the  average  inventory  level 
as  the  simple  average  of  the  initial  cycle  inventory  and  any  residual  inventory.  The  second 
approximation,  the  average  inventory  level  is  simply  plus  the  safety  stock.  The 
assumption  behind  this  approach  is  that  inventory  consists  of  two  separate  parts,  safety  stock 
which  is  always  there,  and  the  regular  inventory  which  is  depleted  through  the  cycle  and 
represented  by 

To  compare  the  two  approximation's  accuracy,  we  have  used  the  crystal  ball 
spreadsheet  simulation  package  to  simplify  the  solvability  of  this  sophisticated  probability 
integral  problem. 

B.  CONCLUSIONS 

As  is  true  with  most  logistical  problems.  Inventory  Officers  in  TALA  are  faced  with 
complicated  demand  forecasting.  Before  advanced  forecasting  methodologies  were  proven 
to  be  effective  they  have  preferred  the  simple  and  traditional  forecasting  approaches  to 
approaches  that  are  puzzling  or  confusing.  Through  this  research,  we  can  make  the  statement 
that  there  is  no  one  specific  forecasting  methodology  that  can  be  applied  universally.  The 
best  way  to  make  demand  forecasting  accurate  and  efficient  is  to  adopt  a  suitable  approach 
that  fits  the  demand  pattern.  In  this  research  we  have  collected  empirical  data  to  explore  the 
various  forecasting  methodology.  The  results  show  that  arithmetic  average  and  moving 
average  are  adequate  for  the  flat  demand  pattern;  however,  more  sophisticated  forecasting 
methods  (EWMA)  are  required  when  demand  data  contains  seasonal  or  trend  components. 
This  result  would  seem  to  require  that  data  must  be  kept  on  each  item  in  the  system  and  the 
method  chosen  for  forecasting  must  fit  the  demand  pattern  for  that  item. 

In  comparing  the  two  holding  cost  area  approximations  with  the  Crystal  Ball 
simulation  package,  the  second  approximation  demonstrated  clear  superiority  and  is 
recommended  for  use  in  modeling  with  appropriate  inventory  policies. 
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C.  RECOMMENDATIONS 


That  TALA  implement  advanced  demand  forecasting  methodologies.  The  advanced 
forecasting  methodologies  discussed  in  Chapter  II  and  III  have  demonstrated  an  improvement 
in  accuracy  over  the  methods  currently  used  by  TALA.  Since  real  demand  pattern  changes 
seasonally  and  due  to  the  effect  of  trends  and  random  factors,  the  current  methods  used  by 
TALA  are  inadequate  and  significant  costs  could  be  saved  if  TALA  were  to  change  their 
forecasting  approach.  Additional  accuracy  could  be  realized  if  holding  costs  and  Ballou 
approximations  of  holding  area  calculations  were  performed  by  TALA. 

Future  research  by  TALA  should  be  performed  since  the  demand  pattern  varies  over 
time  and  it  is  necessary  to  continually  update  the  smoothing  constants  used  in  the  EWMA 
demand  forecasting  method.  At  this  time  updating  smoothing  constants  can  only  be  done  by 
trial  and  error.  It  may  be  possible  for  TALA  to  develop  computer  software  to  calculate 
optimal  smoothing  constants  based  on  their  supply  system  demand  data.  This  would  be  more 
efficient  and  more  accurate  than  hand  calculations  and  is  likely  to  be  more  timely. 
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Demand  Demand 


APPENDIX  A.  FLAT  DEMAND  PATTERN 


»  Actual  Demand  »  Arithmetic  Average  -  TALA's 
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Figure  15.  Flat  Demand  Pattern  -  Arithmetic  Average 


^  Actual  Demand  Arithmetic  Average  -  TALA's 

Last  Period  Forecast  ^  Fypnnentiai  Smoothing,  a=  0.1 
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Figure  16.  Flat  Demand  Pattern  -  Exponential  Smoothing  a=0.1 
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Demand  Demand 


Actual  Demand  ,  Arithmetic  Aserage  -  TALA's 

^  Last  Period  Forecast  Exponential  Smoothing,  a=  0.3 
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Figure  17.  Flat  Demand  Pattern  -  Exponential  Smoothing  a=0.3 

^  Actual  Demand  Arithmetic  Average  -  TALA's 

^  Last  Period  Forecast  ^  Exponential  Smoothing,  a=  0.5 
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Figure  18.  Flat  Demand  Pattern  -  Exponential  Smoothing  a=0.5 
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Demand 


Actual  Demand  _b_  Arithmetic  Awrage  -  TALA's 

_^_Last  Period  Forecast  _h—  Exponential  Smoothing,  a=  0.8 


Figure  19.  Flat  Demand  Pattern  -  Exponential  Smoothing  a=0.8 
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Demand  Demand 


APPENDIX  B.  INCREASING  DEMAND  PATTERN 


,  Actual  Demand  ,  Arithmetic  Average  -  TALA's 


Figure  20.  Trend  Demand  Pattern  -  Arithmetic  Average 


Actual  Demand  ^  3  Month  Moving  Average  -  TALA's 


Figure  21.  Trend  Demand  Pattern  -  Three  Month  Moving  Average 
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Demand  Demand 


..^■■■■Actual  Demand 
^■Arithmetic  Average  -  TALA'S 

^  _  Forecast  Demand  -  Trend  Smoothing  a=0.1,  b=0.1,  T(1)=0 
H  Forecasted  Demand  -  TALA’s  -  3  Month  Moving  Average 


Figure  22.  Trend  Demand  Pattern  -  Trend  Smoothing  a=0.1,  b=0.1,  t(l)=0 


»  Actual  Demand 

Arithmetic  Average  -  TALA's 

-  Forecast  Demand  -  Trend  Smoothing  a=0.5,  b=0.1,  T(1)=0 
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Figure  23.  Trend  Demand  Pattern  •  Trend  Smoothing  a=0.5,  b=0.1,  t(l)=0 
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Demand  Demand 


^  Actual  Demand 
,  Arithmetic  Average  -  TALA’S 

>  Forecast  Demand  -  Trend  Smoothing  a=0.8,  b=0.1,  T{1)=0 
IJZ  Forecasted  Demand  -  TALA’s  -  3  Month  Moving  Average 


Figure  24.  Trend  Demand  Pattern  -  Trend  Smoothing  a=0.8,  b=0.1,  t(l)=0 


Actual  Demand 
Arithmetic  A\«rage  -  TALA's 

^  Forecast  Demand  -  Trend  Smoothing  a=0.8,  b=0.2,  T(1)=0 
w  Forecasted  Demand  -  TALA's  -  3  Month  Moving  Average 


Figure  25.  Trend  Demand  Pattern  -  Trend  Smoothing  a=0.8,  b=0.2,  t(l)=0 
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Demand  Demand 


Actual  Demand 
Arithmetic  Average  -  TALA's 

Forecast  Demand  -  Trend  Smoothing  a=0.8,  b=0.5,  T{1)=0 
-  Forecasted  Demand  -  TALA's  -  3  Month  Moving  A\«rage 


Figure  26.  Trend  Demand  Pattern  -  Trend  Smoothing  a=0.8,  b=0.5,  t(l)=0 


4  .  Actual  Demand 
a  Arithmetic  Average  -  TAIL’S 

^  Forecast  Demand  -  Trend  Smoothing  a=0.8,  b=0.8,  T(1)=0 
H-  Forecasted  Demand  -  TAU\'s  -  3  Month  Movng  Average 


Figure  27.  Trend  Demand  Pattern  -  Trend  Smoothing  a=0.8,  b=0.8,  t(l)=0 
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Demand 


Actual  Demand 

_p_  Arithmetic  Average  -  TALA's 

Forecast  Demand  -  Trend  Smoothing  a=0,5,  b=0.8,  T(1)=0 
Forecasted  Demand  -  TALA’s  -  3  Month  Moving  Average 


Figure  28.  Trend  Demand  Pattern  -  Trend  Smoothing  a=0.5,  b=0.8,  t(l)=0 
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Demand  Demand 


APPENDIX  C.  SEASONAL  DEMAND  PATTERN 


Actual  Demand 
Arithmetic  AMsrage  -  TALA's 
^  Forecast  Demand  -  a=0.1.  c=0.1 
^  -  Forecasted  Demand  -  TALA's  -  3  Month  Moving  Avgrage 


Month 


Figure  29.  Seasonal  Demand  Pattern  a=0.1,  c=0.1 


^  Actual  Demand 

Arithmetic  Average  -  TALA's 
Forecast  Demand  -  a=0.5,  c=0.1 
X  - Forecasted  Demand  -  TALA’s  -  3  Month  Moving  Average 


Figure  30.  Seasonal  Demand  Pattern  a=0.5,  c=0.1 
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Demand  Demand 
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Figure  31.  Seasonal  Demand  Pattern  a=0.8,  c=0.1 


^-Actual  Demand 

Arithmetic  A\erage  -  TALA's 
^  Forecast  Demand  -  a=0.8,  c=0.001 
X  Forecasted  Demand  -  TALA's  -  3  Month  Moving  Average 
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Figure  32.  Seasonal  Demand  Pattern  a=0.8,  c=0.001 


■  Actual  Demand 

Arithmetic  Average  -  TALA's 
^  Forecast  Demand  -  a=0.8,  c=0.1 

Forecasted  Demand  -  TALA’s  -  3  Month  Moving  Average 
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Demand  Demand 


^  Actual  Demand 

Arithmetic  A\«rage  -  TALA's 
^  _  Forecast  Demand  -  a=0.8,  c=0.01 
^  Forecasted  Demand  -  TALA's  -  3  Month  Moving  Average 


Month 


Figure  33.  Seasonal  Demand  Pattern  a=0.8,  c=0.01 


—^Actual  Demand 

Arithmetic  Awrage  -  TALA's 
^  Forecast  Demand  -  a=0.5.  c=0.2 
„  Forecasted  Demand  -  TALA's  -  3  Month  Mowng  Average 


Figure  34.  Seasonal  Demand  Pattern  a=0.5,  c=0.2 
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Demand  Demand 


^  Actual  Demand 

Arithmetic  Average  -  TALA's 

Forecast  Demand  -  a=0.5.  c=0.3 

Forecasted  Demand  -  TALA's  -  3  Month  Moving  Average 


Month 


Figure  35.  Seasonal  Demand  Pattern  a=0.5,  c=0.3 


^  Actual  Demand 
,  Arithmetic  Average  -  TALA’s 
^  _  Forecast  Demand  -  a=0.5,  c=0.5 

Forecasted  Demand  -  TALA's  -  3  Month  Moving  Average 


Figure  36.  Seasonal  Demand  Pattern  a=0.5,  c=0.5 
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Demand 


Actual  Demand 
Arithmetic  Average  -  TALA'S 
Forecast  Demand  -  a=0.5,  c=0.8 
^  Forecasted  Demand  -  TALA'S  -  3  Month  Moving  Average 


Figure  37.  Seasonal  Demand  Pattern  a=0.5,  c=0.8 
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Figure  40.  Combination  of  ROP  &  Q  (ROP=l,  AVELT=0.001) 
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Figure  41.  Combination  of  ROP  and  Q  (ROP=1000,  AVELT=1) 
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Figure  43.  Combination  of  ROP  and  Q  (ROP=1000,  AVELT=0.5) 
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Figure  44.  Combination  of  ROP  and  Q  (ROP=500,  AVELT=0.5) 

0.02 
0 

-0.02 
-0.04 
-0.06 
-0.08 
-0.1 
-0.12 


Figure  45.  Average  Leadtime  Coefficient  of  Variation  vs  Inventor  Area 
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